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Transport of carriers in two dimensional graphene at high electric fields is investigated by combin- 
ing semi-analytical and Monte-Carlo methods. A semi-analytical high-field transport model based 
on the high rate of optical phonon emission provides useful estimates of the saturation currents 
in graphene. For developing a more accurate picture, the non-equilibrium (hot) phonon effect and 
the role of electron-electron scattering were studied using Monte Carlo simulations. Monte Carlo 
simulations indicate that the hot phonon effect plays a dominant role in current saturation, and 
electron-electron scattering strongly thermalizes the hot carrier population in graphene. We also 
find that electron-electron scattering removes negative differential resistance in graphene. Transient 
phenomenon as such as velocity overshoot can be used to speed up graphene-based high speed elec- 
tronic devices by shrinking the channel length below 80nm if electrostatic control can be exercised 
in the absence of a band gap. 

PACS numbers: 81.10.Bk, 72.80.Ey 



I. INTRODUCTION 

Monolayer graphene has attracted great amount of in- 
terest due to its novel electronic properties [1-4]. It 
has perfect two-dimensional (2D) geometry and a lin- 
ear band structure near Dirac points. Since its discov- 
ery, a substantial body of work has been done on charge 
transport properties in graphene [5-9]. Electron mobil- 
ity as high as - 120,000 cm 2 /V-s at - 240 K has been 
measured in suspended graphene [7]. At high field bias, 
epitaxial graphene field effect transistors (FETs) on SiC 
substrates have exhibited high current drives, exceeding 
~ 3 mA//im [10]. The current carrying capability of 
graphene, and the limiting scattering mechanisms have 
been the focus of recent research. 

The earliest studies of high-field current drive sug- 
gested that monolayer graphene should be able to deliver 
currents in the range of ~ 4 mA//im by considering the 
intrinsic optical phonon scattering [11]. Shortly there- 
after, current saturation was experimentally achieved in 
FET structures and the saturation current was observed 
to be proportional to the square root of the channel car- 
rier density [12]. Surface-optical (SO) phonon scattering 
rather than intrinsic optical phonon scattering was sug- 
gested as the dominant scattering mechanism for current 
saturation in graphene on Si02 substrates. The effect of 
hot phonons on limiting the current drive in graphene 
was then considered, but found to be masked by other 
scattering mechanisms for graphene on Si02 substrate 
[13]. Subsequent theoretical and experimental work fur- 
ther confirmed that the limiting mechanism for high-field 
transport in graphene on Si02 substrates is SO phonon 
scattering [14-17]. 

The extrinsic scattering from Si02 could in principle 
be diminished by using different substrates, or by sus- 
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pending graphene. The saturation current in intrinsic 
graphene is still poorly understood. More recently, opti- 
cal phonon temperatures as high as 1000 K was reported 
by Raman measurements in graphene sheets under high 
bias [18], and the hot phonon effect was observed in 
experiments [19]. These findings show the necessity of 
considering hot phonon effect on the intrinsic high field 
transport properties of graphene. Moreover, although 
electron-electron (e-e) scattering degrades the low field 
transport properties of graphene [20], its effect on high 
field transport has not been studied yet. As a result, the 
high- field current carrying capacity of intrinsic graphene, 
considering the hot phonon effect and e-e scattering, is 
not clearly understood. 

In this work, high field carrier transport in intrinsic 
graphene is investigated by Monte-Carlo (MC) simula- 
tions, by including hot phonon effects as well as e-e scat- 
tering. Furthermore, transient high-field transport ef- 
fects, such as velocity overshoot and transit times are 
investigated. In the next section, the various scattering 
mechanisms used in MC simulations are described. In 
the third section, MC simulation results are presented, 
followed by discussions of the results. In the last section, 
we draw conclusions and discuss possible future work. 



II. SCATTERING MECHANISMS IN 
GRAPHENE 

At finite temperature, the intrinsic scattering mecha- 
nisms that affect carrier transport are due to electron- 
phonon interactions. The phonon modes that need to be 
accounted for are acoustic and optical. The scattering 
rates are calculated using Fermi's Golden rule [21]. The 
acoustic phonon scattering rate from state \k) to state 
\k') is given by 
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where D ac ~ 16 eV is the deformation potential of acous- 
tic phonon scattering, ks is the Boltzmann constant, T 
is the temperature, Sk = hvp\k\ is the kinetic energy 
of the carrier and k is the wavevector, h is the reduced 
Planck constant, vp = 10 8 cm/s is the Fermi velocity 
in graphene, a m ~ 7.6 x 10 -8 g/cm 2 is the 2D mass 
density of graphene, and v p = 20 km/s is the acous- 
tic phonon (sound) velocity in graphene. Okk' is the an- 
gle of electron velocity between initial state \k) and final 
state \k'). Acoustic phonon scattering in graphene is not 
isotropic, since the acoustic phonon perturbation is long- 
range, and the wave functions of carriers in graphene 
have the 'pseudospin' symmetry. The angle dependent 
factor for acoustic phonon scattering is l + cos^fc 7 ? which 
prevents backscattering. The acoustic phonon scattering 
rate for electron in state \k) is calculated by summing 
over all possible finals in k space and it is given by 
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where acoustic phonon scattering is treated as quasi- 
elastic, and both emission and absorption are considered. 

The optical phonon scattering rate from state \k) to 
state \k') is given by 
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where D Q is the mode-specific optical deformation poten- 
tial of graphene, uj is the optical phonon frequency and 
n op = 1/ '( e huj o/k B T _ ^ g foe optical phonon number. 
Similarly, the optical phonon scattering rate for state \k) 
is the summation of all scatterings to final states, and it 
is given by 
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where ± is for phonon absorption and emission respec- 
tively. In graphene, the zone edge (ZE) transverse opti- 
cal (TO) mode has the strongest coupling with electrons 
[22] and the energy of this mode is around hu = 160 
meV. The 'optical' deformation potential of the mode 
is D Q = 25.6 eV/A [23]. The scattering of zone-edge 
mode is non-isotropic and has the angle dependence fac- 
tor 1 — cos 6kk' [24]. Forward scattering (0 = 0) is pro- 
hibited and back scattering (6 = tt) is preferred due to 
this factor; the situation is opposite of what happens in 
acoustic phonon scattering. In our simulation, we also 
include the zone center (ZC) optical phonon. The ZC op- 
tical phonon has a higher energy, Hlj = 196 meV. The 
deformation potential for the ZC phonon is D Q = 14.1 
eV/A[25]. The zone center optical phonon scattering is 
isotropic [24]. 

The effect of e-e scattering in graphene with linear 
£ ~ k dispersion is different from traditional semicon- 
ductors that have parabolic band structure (£ ~ k 2 ) [20]. 
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FIG. 1: Scattering rates versus energy in 2D graphene: 
T=300 K, n e = 7.7 x 10 12 cm" 2 . The dashed line is the 
sum scattering of acoustic phonon and optical phonons. 



Energy and momentum conservation requirements mod- 
ify the carrier distribution and center of mass motion in 
a way that the low-field transport properties (such as 
carrier mobility) is strongly damped, especially at low 
carrier concentrations. The e-e scattering rate of elec- 
tron at state \k) with another electron at \ko) is given 

by 
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where k' and k' are the finals states. The Coulomb in- 
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where wave vector q = \k — k'\ and e(q) is the dielectric 
function of 2D graphene. O^y and k2k > 2 are the scat- 
tering angles of the two electrons, which has to be taken 
into account since Coulomb scattering is also long range 
interaction. Similarly, wave vector q' = \k — k' \ is for the 
scattering process of exchange of two electrons. The two 
delta function shows the conservation of momentum and 
energy in the scattering. The electron scattering rate is 
calculated by summation of all ko state and also all the 
final states (k',k' Q ), and it is given by 
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where the occupation of state fco is the Fermi-Dirac dis- 
tribution and the final states are assumed to be empty. 
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The e-e scattering rate is incorporated into the MC 
simulation in this work following the procedure of ref. [20] 
with one difference: the screening of Coulombic forces is 
treated as temperature-dependent [26]. 

The scattering rates at 300 K are plotted in Fig.l. 
In order to capture high-field effects, the scattering rate 
from to 1 eV are calculated for all the scattering mecha- 
nisms in intrinsic graphene. The intrinsic optical phonon 
scattering rate increases linearly with energy at states 
higher than huo • Acoustic phonon scattering rate is one 
order of magnitude smaller than the optical phonon scat- 
tering at high energies. Electrons are accelerated by elec- 
tric field and populate high-energy states. Therefore, op- 
tical phonon scattering is the dominant energy relaxation 
channel for hot electrons rather than acoustic phonon. 
Since the ZE optical phonon has lower energy and higher 
scattering rate, it is the dominant optical mode. 

The e-e scattering rate at 300 K is higher than optical 
phonon scattering. This scattering thermalizes the car- 
rier distribution, driving it towards a Fermi Dirac distri- 
bution. The net momentum does not change during e-e 
collisions, but the average velocity of electrons changes 
due to the linear band structure of 2D graphene. As a 
consequence, e-e scattering has an impact on the current. 
The e-e scattering rate increases with carrier temperature 
so that it become more important under high field. 

From Fig.l, the phonon scattering rate at high energies 
exceeds 10 ps _1 . Optical phonons emitted by electrons 
have to relax to other phonon modes by anharmonic in- 
teraction before the heat propagates into the substrate 
or contacts. The characteristic lifetime of optical phonon 
decay into acoustic phonon modes is ~1 to 5ps in carbon- 
based sp 2 crystals [27]. Since the optical phonon lifetime 
is much longer than the phonon generation time (~ 0.1 
ps), the optical phonon population in graphene is out 
of equilibrium. The fast generation rate and slow decay 
rate of optical phonons result in a hot phonon effect in 
graphene. 

Before presenting the numerical results of the ensem- 
ble MC calculation, a heuristic picture of the high field 
transport properties in graphene is developed in the next 
section. A simple model considering the effect of carrier 
degeneracy provides analytical estimation of the satura- 
tion currents for degenerate electron systems. 



III. HIGH FIELD TRANSPORT MODEL IN 
GRAPHENE 

The electron distribution function in the k- space in the 
steady state at high electric fields can be approximated 
since the optical phonon emission is much faster than 
competing energy relaxation mechanisms. The steady 
state electron distribution function /(k) under this as- 
sumption is shown in Fig. 2 (a). The electrons are dis- 
tributed between two energy contour circles, the high 
energy circle Eh = hvpkh and the low energy circle 
£i = hvpki. All electrons move toward k x direction un- 
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from which the dimensions of the distribution function 
(kijkh) can be calculated at any carrier concentration. 
The current flowing in the direction of the electric field 
corresponding to the distribution function in Fig. 2 (a) is 
then found by summing the x— directed group velocities 
of all states inside the shaded region in the k space; the 
value is found to be 
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FIG. 2: Current saturation model in 2D graphene: (a) The 
distribution function of streaming model at steady state. 
f(k) = l at colored region, otherwise, f(k)=0. (b) The satu- 
ration velocity versus carrier concentration. 



der the applied electric field. Once electrons approach the 
high energy circle on the right, they are scattered back 
to low energy circle by emitting optical phonon instan- 
taneously. Therefore, the energy difference between the 
two circles is the optical phonon energy, £ ^ — E\ = hw 01 
as indicated in the figure. The state occupation inside 
the low energy circle is full so that optical phonon emis- 
sion by electrons with energy smaller than £h are prohib- 
ited by the Pauli exclusion principle. As optical phonon 
emission continues and scatters the electrons from high 
energy states back to the low energy states, distribution 
function /(k) is squeezed and elongated along the direc- 
tion of the electric field. The steady distribution function 
is called 'streaming' distribution. This streaming model 
captures the effect of carrier degeneracy, which has not 
been considered in high field-transport in traditional non- 
degenerate semiconductors. 

The carrier concentration is proportional to the area of 
the distribution function S in k space, n = g s g v S/(27r) 2 , 
where g s and g v are spin and valley degeneracies in 
graphene. The area of the distribution function is given 

by 

h 
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where £\ is the electron energy at low energy circle. This 
streaming model shows that the saturation current is pro- 
portional to the optical phonon frequency and radius of 
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low energy circle. In Fig. 2(b), the ensemble saturation 
velocity v sat = J sa t/en is plotted against the electron 
density. The saturation velocity varies with carrier con- 
centration and is not a constant. It approaches the Fermi 
velocity in graphene vp = 10 8 cm/s as n —> 0, and falls to 
~2 x 10 T cm/s at n ~ 10 13 /cm 2 . The dotted line shows 
the closely related square root relation v sat = cj /2y / n7r 
proposed in ref.[ll]. 

Traditionally, the saturation velocity in non- 
degenerate semiconductors has been considered a 
constant, e.g. v sat = 10 7 cm/s in Si at room tempera- 
ture [28]. The reason for the carrier-density dependent 
saturation velocity in graphene is the degeneracy of 
carriers and the Pauli exclusion principle. Close to 
inverse square-root dependence of the effective ensemble 
saturation velocity on the carrier density has been mea- 
sured in graphene on Si02 substrate [12, 17]. The model 
presented here captures the essence of the high-field 
transport and is in agreement the experimental results, 
indicating that the carrier degeneracy plays a major role 
on high field transport in graphene. However, the model 
here can stand on a firmer footing by using a numerical 
Monte- Carlo approach to find high field transport 
properties, which is the topic of the next section. 



IV. MONTE CARLO SIMULATION: 
IMPLEMENTATION AND RESULTS 

In the MC implementation of the high-field transport 
problem, we have considered electron energies ranging 
from to 10 hu in the graphene conduction band. The 
k— space is divided into around 10 5 cells. In simulation, 
'superelectrons' move in k space according to Newton's 
law F = hdk/dt [21]. The number of 'superelectrons' 
varies from 1000 to be more than 10 4 , depending on the 
carrier density. The current is calculated by summation 
of the group velocities over all the 'superelectrons'. 

The screening of the e-e Coulombic interaction is de- 
pendent on the carrier temperature. In the MC simula- 
tion, the electron temperature T e is updated every time 
slot (~0.5 fs). The temperature was calculated from the 
thermal energy of electron gas £ t h = — 
where hi is wavevector of the i'th electron, and the av- 
erage wavevector of all electrons is k. The electron tem- 
perature was calculated by assuming a Fermi-Dirac dis- 
tribution, whereby the thermal energy of the electron gas 
at temperature T e is given by 
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where A is the area of graphene sheet and r] = £f/kBT e . 
In simulation, we assume the carrier density in graphene 
conduction band is fixed by ignoring generation pro- 
cesses. As a result, Fermi level changes with the tem- 
perature, and both of them can be calculated once we 
know the thermal energy and carrier concentration. 
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FIG. 3: (a) The current density versus electric field at two 
carrier concentration nl=7.5 x 10 12 cm -2 and n2=3.0 x 10 12 
cm -2 . Open symbols: no hot phonon; Solid symbols: hot 
phonon (t p h — 1 ps). Dot-dash lines consider hot phonon ef- 
fect but not e-e scatterings, (b) Saturation current versus car- 
rier density. Solid black curve: the current saturation results 
of streaming model; Solid square: intrinsic graphene without 
hot phonon effect, T — 50 kV/cm; Open square: with hot 
phonon effect (r p h=l ps), T — 50 kV/cm. All simulations 
are at 300 K. 



In order to capture the effect of hot phonons, the 
phonon q-space is discretized in the same fashion as the 
electron k- space, and the hot-phonon effect is treated as 
described in [29]. The optical phonon distribution func- 
tion is updated based on the phonon wavevectors dur- 
ing the emission and absorption events. Optical phonon 
decay into other phonon modes is governed by the rate 
equation of phonon occupation function n q (t) - 



n q (t + St) = n q (t) 



n q (t) - n q o 



St, 
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where St is the time slot in MC simulation, n q o is the equi- 
librium Bose-Einstein distribution, and r p h is the phonon 
lifetime. The time slot St is around 0.5 fs, correspond- 
ing to a scattering rate 2xl0 3 ps -1 , which is higher than 
the summation of all scattering rates for electrons. This 
time slot was updated according to the phonon popu- 
lation in simulation, since the maximum scattering rate 
varies with the phonon population. 

The degeneracy of carriers is explicitly dealt with by 
the 'rejection technique' [30]. The simulation is continued 
till a steady state distribution function is obtained. This 
procedure yields the distribution functions of electrons 
f(k x ,k y ,t) and phonons n(q x ,q y ,t) at every time step, 
which is then used to find ensemble velocities, currents, 
as well as transient effects such as velocity overshoot as 
a function of the carrier density and the applied electric 
field. 

The steady-state currents calculated by the MC sim- 
ulation for various strengths of electric fields are shown 
in figure 3(a) for two carrier densities. The curves of dif- 
ferent carrier concentrations show similar behavior. The 
currents increase with the field and saturate at a field 
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T ~ 50 kV/cm. The degradation of the saturation cur- 
rent due to the hot-phonon effect is clearly seen by com- 
paring the open and solid symbol lines, which is indi- 
cated by the arrows. Hot phonon effect is more signif- 
icant in degradation of the saturation current at higher 
carrier concentrations. In order to highlight the effect 
of e-e scattering at high field, the transport without e- 
e scattering is also simulated. The dot-dash lines are 
the currents including hot phonon effect but without e- 
e scatterings. Without e-e scattering, the currents drop 
at high fields, which is the indication of negative differ- 
ential resistance (NDR). NDR has been proposed and 
observed in simulations for carbon naotubes (CNTs) and 
graphene [31-33], due to the linear dispersion of band 
structure. In graphene, the electrons at high energy state 
move at the same speed with the low energy state, which 
is the constant Fermi velocity in graphene. At high field, 
more electrons move into high energy states. But the 
increasing of velocity is not as rapidly as conventional 
parabolic materials. As more and more electrons occupy 
high energy states the backscattering is increased, which 
degrades the average velocity of electron gas and leads 
to NDR. However, this effect disappears (or weakened) 
if e-e scattering is considered. The e-e scatterings ex- 
change electrons' momentum and energy, diving electrons 
into Fermi-Dirac distribution. As a result, high energy 
electrons exchange their momentum and energy with low 
energy electrons before substantial backscatterings hap- 
pen. With e-e scatterings, the hot electrons generation 
and backscattering are weakened and the ensemble ve- 
locity of carriers increases with the applied electric field. 
The current-field curves from the MC simulation there- 
fore show that a) the hot phonon effect acts as a strong 
current limiting mechanism at high carrier densities and 
b) e-e scattering gets rid of NDR effect in graphene. 

In Fig. 3(b), the saturation current is plotted against 
the carrier density. The solid black curve is the analyt- 
ical result from streaming model (equation 9) for com- 
parison. The solid squares are the intrinsic saturation 
currents without the hot phonon effect. The saturation 
current is slightly lower than the analytical result due to 
hot phonon degradation. Saturation current density of 
~3 A/mm has been measured in epitaxial graphene on 
SiC substrate (n ~ 10 13 cm -2 ) [10], which is consistent 
with the MC result. The current-carrier concentration 
curves are not linear, indicating the ensemble saturation 
velocity is not constant, but decreases with increasing 
carrier density. 

The hot phonon effect is better illustrated by the distri- 
bution functions of electrons in k- space, and phonons in 
q-space. The electron distribution functions are plotted 
in Fig. 4(a) along the k x axis which is along the direc- 
tion of the electric field. The solid line is the electron 
distribution at steady state, considering hot phonon and 
e-e scattering. The peak of the distribution function is 
close to 1, showing the high degeneracy of carriers in 
graphene. For comparison, we turned off the e-e scatter- 
ing or hot phonon effect in the simulation. The distribu- 
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FIG. 4: (a) The electron distribution function along k x with 
k y =0. Solid line considers hot phonon and e-e scattering. 
(Dash line: turning off hot phonon effect and dash-dot line: 
turning off e-e scattering in simulation) (b) The phonon 
number distribution function in q space. Carrier density: 
n=7.5xl0 12 cm -2 , electric field T — 50 kV/cm, phonon life- 
time r p h=l ps, environment temperature T =300 K. 



tion function without e-e scattering has a lower peak and 
is broader, shown by the dot-dash curve. Compared to 
the solid distribution curve, the distribution without e-e 
scattering has more electrons at high energy states, but 
also more backscattered electrons. This highlights the 
role of e-e scattering in maintaining the shape of carrier 
distribution close to Fermi-Dirac like (the solid line) and 
thus preventing the NDR effect. Similarly, the effect of 
hot phonons becomes clear when we compare the distri- 
bution functions shown by the solid and dashed curve. 
The carrier distribution with hot phonons is shifted less 
along the field direction, so the saturation current is low- 
ered by the hot-phonon effect. 

The effect of hot phonons on the transport can also be 
seen from the phonon distribution in q-space. Fig. 4(b) 
shows the 2D phonon distribution functions n(q xi q y ). 
The hot phonon distribution has a peak at positive 
q x « buoo/vF-, which corresponds to the fast optical 
phonon emission by hot electrons. The phonon occu- 
pation number peak exceeds n p h ~ 5 which corresponds 
to a phonon temperature of T op ~ 10 4 K. 

Besides the steady-state saturation currents, the MC 
simulation procedure also allows us to study transient 
phenomena in high- field carrier transport. Prior to ap- 
proaching steady state, carriers moving short distances 
in traditional semiconductors are known to undergo tran- 
sient velocity overshoot, attaining higher velocities than 
what can be attained in the steady state. The hot elec- 
trons accelerated by high field keep moving in k space 
without significant phonon emission during timescales 
shorter than the mean phonon emission time. As a re- 
sult, a velocity peak is achieved in a short period. In high 
speed electronic devices, the channel length is scaled to 
tens of nanometers. The carriers travel from source to 
drain contacts within ps timescales, and velocity over- 
shoot can occur. This transient velocity overshoot can be 
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utilized to lower the transit time of carriers in graphene 
devices, thus resulting in an enhanced speed. Fig. 5 (a) 
shows the ensemble transient velocity response to electric 
fields. Velocity overshoot occurs when the electric field 
is higher than 5 kV/cm. The higher the field, the higher 
is the overshoot velocity peak; it approaches ~ 5 x 10 T 
cm/s at a field of 50 kV/cm for the chosen carrier density 
(n ~ 7.5 x 10 12 /cm 2 ). After reaching the peak overshoot 
velocity, the ensemble drift velocity decreases gradually 
by dissipating excess energy into the graphene crystal in 
the form of optical phonon vibrations. The long tail of 
decreasing velocity is due to the hot phonon effect - op- 
tical phonons are slow in decaying into acoustic modes 
which heat up the lattice. 

Velocity overshoot occurs within a time window of 
t os ~ 0.2 ps. If the channel length (L) is short enough 
(L < L os = J tos Vd(t)dt), the hot carrier energy is dissi- 
pated at the drain contact, and the hot phonon effect can 
be prevented in a short channel device. Thus, the dis- 
tance L os = J^ os Vd(t)dt that carriers travel in the veloc- 
ity overshoot period is crucial for the design of high-speed 
devices. The transit times t tr given by L = J^ tr Vd(t)dt 
is calculated from the result in Fig5(a) and shown for 
various lengths in 5(b). For short channel lengths, the 
transit time is sensitive to the magnitude of the electric 
field, since velocity overshoot increases with electric field. 
For long channel lengths, a major portion of the transit 
time is spent traveling at the steady state drift velocity, 
and the overshoot velocity makes a minor contribution. 
In the MC simulation, the transit time is not found to be 
particularly sensitive to the carrier density, which bodes 
well for devices requiring high currents and high speeds. 
In a L = 80 nm long channel device, the transit time 
is found to be t tr ~ 0.2 ps, which is exactly the time 
window for velocity overshoot. Thus, graphene channels 
L < 80 nm can fully take advantage of velocity overshoot 
for enhancing their speed. 

The peak overshoot velocity, and overshoot time win- 
dow (t s) is determined by strength of electron-optical 
phonon interaction, the optical phonon energy, as well as 
the low-field mobility. The optical phonon energy and 
carrier-phonon interaction are intrinsic properties of ma- 
terial, but the peak overshoot velocity peak can be im- 
proved by increasing the cleanliness of the substrate or 
by using gate dielectrics with less impurities. 

V. DISCUSSION AND CONCLUSION 

We point out some limitations of this work, and suggest 
topics that need investigation in the future. In this work 
we have ignored the band to band processes, e.g. inter- 
band (Zener) tunneling and e-e scattering induced carrier 
generation. The carrier concentration does not change at 
high electric field under the assumption. However, this 
is not the case in graphene if the Fermi level is close to 
the Dirac point. Since the band gap is zero, tunneling 
and e-e scatterings generate electrons and holes under 
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FIG. 5: (a) Electron drift velocity evolution with time under 
field (b) Carrier transit time versus channel length. Optical 
phonon lifetime r p ^=l ps, environment temperature T=300 
K, carrier density n — 7.5 x 10 12 cm -2 . 



high fields at very low carrier concentrations. In exper- 
imental measurements, the current increases and never 
saturates when the Fermi level is close to the Dirac point 
in intrinsic graphene (Sf ~ 0) [13]. Field-induced car- 
rier generation is also the reason of the low on/off ratio 
(~ 10) in graphene FETs. 

In experiments, intrinsic properties of graphene is dif- 
ficult to achieve. SO phonon scattering limits the high 
field transport for graphene in close proximity to sub- 
strates with low energy phonon modes, like Si02- Sus- 
pended graphene is free from environmental scattering 
effects, but the carrier density is not easy to modulate. 
One possible way to measure the intrinsic performance is 
to use SiC substrate graphene which has weak SO phonon 
coupling. 

We assumed the lattice temperature to be at 300 K in 
our simulations, although the optical phonon tempera- 
ture can be much higher than environment temperature 
300 K. This assumption could overestimate the satura- 
tion current, since the lattice temperature in graphene 
increases due to the limited thermal conductance of the 
substrate and contacts. The hot phonon effect is dom- 
inant, so to a certain extent ignoring the rise of lattice 
temperature is justified. 

The phonon lifetime has been measured for graphene 
on SiC substrate and it was found to be around 2.5 ps 
[34]. The hot phonon effect can be diminished by in- 
troducing isotopic disorder into graphene [35]. Isotopic 
disorder introduces phonon modes that are localized in 
real space, and therefore spread in q- space, and more 
phonon modes are thus involved in the cooling of car- 
riers. The saturation currents can be improved by this 
disorder engineered lifetime of optical phonons, with di- 
rect consequences of higher saturation currents. 

In conclusion, we investigated the high field transport 
in intrinsic monolayer graphene. MC simulation results 
show that both hot-phonon effect and e-e scattering have 
strong impacts on the saturation current. Transient ve- 
locity overshoot was studied and the short carrier transit 
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